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ABSTRACT 


The  electromagnetic  surface  currents  induced  by  a  magnetic 
dipole  on  infinite  perfectly  conducting  surfaces  in  the  presence  of 
a  sea-air  interface,  as  well  as  in  a  homogeneous  ocean  environ¬ 
ment,  are  obtained.  The  method  used  is  that  of  images  where  the 
approximation  of  a  -1  reflection  coefficient  at  the  sea-air  interface 
may  be  justified  and  used.  Computer  programs  and  results  which 
yield  a  clear  insight  into  the  sphere  of  influence  of  the  magnetic 
dipole  in  such  an  environment  when  displayed  graphically  are  given. 
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ELECTROMAGNETIC  SURFACE  CURRENTS  INDUCED  BY  A  MAGNETIC 
DIPOLE  SOURCE  ON  INFINITE  PERFECTLY  CONDUCTING 
SURFACES  IN  THE  PRESENCE  OF  AN  AIR-OCEAN  INTERFACE 


INTRODUCTION 

In  the  first  part  of  this  report,  the  electromagnetic  surface  current 
induced  on  a  perfectly  conducting  infinite  flat  plate  in  a  conducting  medium  for 
a  magnetic  dipole  source  is  obtained.  These  results  will  allow  the  introduction 
of  the  effect  of  the  sea-air  interface  in  the  second  part  of  the  report.  The 
analysis  employs  the  method  of  images. 

The  results  will  serve  to  illustrate  the  extent  of  the  dipole  influence  when 
immersed  in  a  conducting  fluid,  and  thus,  evaluate  the  effectiveness  of  the 
magnetic  dipole  as  a  detector  in  the  ocean. 


THEORY 

INDUCED  CURRENT  ON  INFINITE  FLAT  PLATES  IN  A 
LOSSY  ENVIRONMENT:  MAGNETIC  DIPOLE  EXCITATION 

The  electric  field  due  to  a  magnetic  dipole  source  is  given  by 

EVEo  =  “j[  (M  2  +  j(kr)  *]  e  ^  sin  9  a^.  (1) 

Thus  the  total  field  at  P  due  to  source  and  image  is 

— jkr 

VE0  =  ^{[K)"2 + W1] e  1  sin  ei 

™  jkr 

-  [K)"2  +  W1] e  2  sin  6  2} a*  (2) 

where 

2  2 

E  =  wuNIdAk  /4 tr=  uf, ik  |ml/4ir. 


(3) 
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NIdA  is  the  magnitude  of  the  magnetic  dipole  moment;  N  is  the  number  of  turns 
on  the  loop.  Also 

k  =  (1  -  j)/4  (4) 

1/2 

4  =  (2  /who)'  .  (5) 

Figure  1  illustrates  the  geometry  considered,  while  figure  2  shows  the  solution 
using  the  method  of  images.  Appendix  A  represents  the  computer  program 
used. 


The  analysis  is  expanded  to  include  the  case  where  the  dipole  orientation 
with  respect  to  the  normal  to  the  conducting  plate  is  arbitrary  and  the  approach 
to  the  plate  may  be  at  any  angle  0p.  The  following  results  pertain  to  the  two 
cases:  a  magnetic  dipole  source  normal  to  the  path,  and  a  magnetic  dipole 
source  in  the  direction  of  the  path.  For  both  cases,  the  current  pattern  induced 
on  the  conducting  plate  and  the  signal  induced  in  the  detector  coil  will  be  calcu¬ 
lated.  Figure  3  illustrates  the  generalized  geometry,  while  figure  4  illustrates 
the  image  theory  solution  for  the  second  case  when  0p  =  */2. 

A  general  analysis  which  can  handle  all  possible  cases  will  be  used,  where 
the  dipole  source  is  given  arbitrary  orientation  in  space  and  the  resultant  field 
at  any  point  in  the  same  half  space  is  obtained.  A  particular  calculation  will 
thus  need  to  have  the  dipole  orientation  and  location,  as  well  as  the  detector 
location  and  orientation  with  respect  to  the  dipole  source,  specified.  A  path 
for  the  dipole  source-detector  combination  will  be  specified  through  controlling 
successive  dipole  source -detector  locations  along  a  linear  path  or  along  any 
general  curve.  In  this  study  primary  interest  lies  in  obtaining  the  induced  sur¬ 
face  currents  on  the  conducting  plane  at  y  =  -d. 

Since  the  dipole  field  components  are  given  by 

Hr/H0  =  2  [(kr)  3  +  j(krf 2]  e’^  cos  B  (6) 

q  _2  “1  -ikr 

VHo  =  [(kr)  +  j(kr)  -  (kr)  ]  e  J  sin  B  (7) 

E^/Hq  =  -j(w#iA)  [(kr)  2  +  j(kr)  1]  e  ^  sin  6  (8) 
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Figure  3.  Generalized  Geometry 
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where 

H  =  NIdAk3/4ir  (9) 

o 

and 

NIdA  =  magnitude  of  magnetic  dipole  moments,  (10) 

in  terms  of  a  local  spherical  coordinate  system,  this  system  must  be  related  to 
the  rectangular  coordinate  system  which  will  describe  the  various  bounding  sur¬ 
faces.  The  magnetic  dipole  source  will,  therefore,  be  taken  to  be  oriented  with 
angles  0S,4>S  with  respect  to  the  rectangular  system  x,y,z,  with  the  z  axis 
coincident  with  0S  =  0  and  the  x  axis  coincident  with  <t>s  =  0.  The  perfectly 
conducting  plane  is  taken  to  be  the  plane  y  =  -d.  Figure  5  illustrates  this  gen¬ 
eral  situation. 

Consider  the  point  (x'.y'.z')  which,  in  the  local  dipole  spherical  coordi¬ 
nate  system,  is  the  point  (r,  0,4>);  then  the  relationship  between  the  two  coor¬ 
dinate  systems  is  given  by 

_  2  2  2-,  1/2 

r  =  0'  +y'  +Z'  ]  /  (11) 

sin  8  =  (x'2  +  y'2)1//2/r  (12) 

cos  8  =  z’/r  (13) 

where 

x'  =  x  cos  <t>  +  y  sin  <t>  (14) 

s  a 

y'  =  f-x  sin  <t>  +  y  cos  <t>  )  cos  8  -  z  sin  8  (15) 

'  s  a'  s  s 

z’  =  (-x  sin  <t>  +y  cos  <fi  )  sin  8  +  z  cos  8  .  (16) 

Thus,  given  the  location  (x,y,z)  at  which  the  dipole  field  needs  to  be  evaluated 
and  the  angles  (8a,<t>a)  giving  the  direction  of  the  dipoles,  the  field  components 
can  be  calculated  using  (14)  through  (16)  in  sequence  with  (11)  through  (13). 
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In  this  case,  because  only  infinite  flat  conducting  planes  are  considered, 
the  induced  surface  current  is  given  by 

J  =  2n  x  H  |  (17 


which  is  obtained  from  either  physical  optics  or  the  theory  of  images,  as  these 
two  approaches  are  fully  equivalent  in  this  situation.  This  can  be  rewritten  as 


J  =  2a  x 

~  ~y 


HI  =  2a  x  (H  a  +  H  aj 

“Ud  -y 


and  the  relationship  between  the  unit  vectors  in  the  local  dipole  spherical  sys¬ 
tem  and  the  rectangular  coordinate  system  must  be  established.  This  is 

a  =  a  fcos  4>  sin  8  cos  <b  -  sin  <p  cos  8  sin  8  sin  <t>  -  sin  <t>  sin  8  cos  (T| 

~r  -x L  s  ss  ssJ 

+  a  fsintf)  sin  8  cos  <t>  +  cos  <t>  cos  8  sin  8  sin  <t>  +  sin  8  cos  0  cos 

^yL  s  s  s  s  s  J 

+  a  T-sin  8  sin  8  sin  <t>  +  cos  8  cos  (19) 

~z L  s  s  J 

a.  =  a  fcos  4>  cos  8  cos  <t>  -  sin  <t>  cos  8  cos  8  sin  <t>  +  sin  8  sin  <t>  sin  fll 

~8  ~xL  s  s  s  s  s  J 

+  a  Tsin  <i>  cos  8  cos  +  cos  8  cos  <b  cos  8  sin  <fi-  sin  8  cos  <t>  sin 

~yL  s  s  s  s  s  J 

+  a  [-sin  8  cos  8  sin  <t>-  cos  8  sin  9]  .  (20) 

z  s  s 

The  two  components  of  the  surface  current  density  then  take  the  form 

J  =  <2H  (cos  6  cos  8  -  sin  8  sin  ^  sin  8  ) 
x  l  r\  s  s' 

-2^  (sin  8  cos  8  +  cos  8  sin  4>  sin  8  )j  (21) 

8  y=-d 
and 

-J  =  <2H  sin  8  (cos  $  cos  $  -  «in  0  cos  8  sin  -  cos  8  sin  ^  sin  0 
z"rL's  s  S'  s  sj 

+  2H^  fcos  8  (cos  p  cos  0  -  sin  0  cos  8  sin  <$)  +  sin  8  sin  <*>g  sin  8 1 11  (22) 
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where 


,  ,2  ,2\l/2 

sin  <t>  =  y'/(x'  +  y'  ) 

(23) 

,  ,/  ,2  2\l/2 

cos  <t>  =  x’/(x'  +  y  j 

(24) 

supplement  the  earlier  results  (11)  through  (13). 

Figure  6  gives  sample  calculations  for  a  dipole  oriented  normal  to  a  con¬ 
ducting  plane.  Because  of  the  rotational  symmetry,  the  current  density  is 
plotted  in  terms  of  the  cylindrical  component  rather  than  the  rectangular 
components  given  in  (21)  and  (22).  The  plot  is  normalized  with  respect  to  a 
reference  surface  current  magnitude  as  shown  on  the  ordinate  axis.  The 
parameter  for  the  various  curves  is  the  distance  of  the  dipole  from  the  plane  in 
tenths  of  one  skin  depth.  It  is  to  be  observed  from  figure  6  that,  even  with 
increasing  separation  in  the  ocean  between  source  and  conducting  plate,  where 
the  electrical  parameters  have  been  taken  to  be  6r  =  80,  a  =  4,  the  region  of 
influence  does  not  extend  significantly  beyond  a  region  given  by  a  circular  area 
of  one  skin  depth  in  radius  centered  at  the  point  on  the  plane  nearest  the  dipole . 


MODIFICATION  OF  THE  INDUCED  CURRENT  DISTRIBUTION 
ON  VERTICAL  IN  FINITE  CONDUCTING  PLATE  DUE 
TO  THE  PRESENCE  OF  AIR -OCEAN  INTERFACE: 

MAGNETIC  DIPOLE  EXCITATION 

Approximate  boundary  conditions  will  be  used  in  the  sense  that  the  con¬ 
ductivity  of  the  conducting  plate  will  be  assumed  infinite  (n  x  E  =  0,  n  •  H  =  0) 
while  at  the  air-ocean  interface  the  dielectric  discontinuity  is  assumed  suffici¬ 
ently  large  so  that  the  approximation,  nxH  =  0,  n-E=0  can  be  used.  In 
terms  of  reflection  coefficients  this  means  that  at  the  plate  p  =  +1  while  at  the 
air-ocean  interface  p  =  -1.  Under  this  approximation,  the  image  theory  can 
again  be  used  to  yield  a  solution  in  a  direct  manner .  It  is  to  be  expected  that 
the  results  obtained  should  be  good  unless  the  dipole  source  is  very  close  to  the 
air-ocean  interface.  The  approximation  at  the  metal  plate  will  be  excellent  for 
any  metal  in  common  usage,  provided  the  thickness  of  the  metal  is  of  the  order 
of  the  characteristic  skin  depth  in  the  metal.  At  the  air-ocean  interface,  where 
at  the  frequencies  considered  here  the  discontinuity  in  the  relative  dielectric 
constant  or  relative  permittivity  is  about  80,  any  ray  incident  on  that  interface 
from  below  will  suffer  total  internal  reflection  unless  the  angle  that  the  ray 
makes  with  the  normal  to  the  surface  is  less  than  approximately  6*. 
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ure  6.  Surface  Current  Density  Induced  on  the  Plane  Normalize 
ch  Respect  to  the  Surface  Current  at  1  Skin  Depth  from  Nearest 
Point  for  1  Skin  Depth  Stand-off 
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Within  this  small  cone,  the  amount  of  energy  transmitted  into  the  air  region 
will  be  extremely  small.  This  is  due  to  both  the  small  area  such  a  cone  projects 
upon  the  surface  for  sources  within  a  few  skin  depths  from  the  surface  and  the 
dielectric  discontinuity  which  yields  a  transmission  coefficient  of  approximately 
0.2,  i.e. ,  a  power  transmission  coefficient  of  0.04.  The  above  approximations 
can,  therefore,  be  used  with  considerable  confidence. 

Figure  7  shows  the  image  theory  solution  for  a  dipole  oriented  vertically. 
The  coordinate  axes  of  the  cartesian  system  are  oriented  such  that  the  x-y 
plane  lies  on  the  air-ocean  interface  with  the  origin  at  the  conjunction  of  the 
air-ocean  and  metal-ocean  interfaces,  the  y  axis  being  placed  normal  to  the 
metal  surface. 

The  components  of  the  total  magnetic  field  intensity  at  any  point  in  the 
ocean  quadrant  are  given  by 


4  — 

VHo  ■  S[»K)'S  *  3iK)'2  -  hr1] e  ‘ 

4  -jkp 

VHo  ■  X^h)-3  *  3ih)'2  -  h)'1] e 

i=l 

Hz/Ho  =  Z  {[2H)'3  +  2i(kri)"2]  cos2  ei 


sin  6.  cos  8.  cos  <t>. 
l  l  l 


sin  8.  cos  8.  sin  <t>. 
l  l  l 


(25) 


(26) 


—  ijcj* 

-  [(kr.)“3  +  j(kri)'2  -  (kr.)-1]  sin2  tfje  1 

(27) 

where,  on  the  surface  y 

=  0, 

r4=ri 

*4=91  *4“  2t  " 

(28) 

r3  =  r2 

93=92  «3=2’r-^2=-^2  =  ^r 

(29) 
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On  this  surface,  the  induced  surface  current  is  given  by 

J  =  n  x  H 


Therefore, 


=  a  x  ("h  a  -  Hal 
-y  [  x~x  z~zj 


=  H  a  -Ha  . 
zrx  x~z 


J  /H  =  H  /H 
X  o  z  o 

J  /H  =  -H  /H 
Z  O  X  o 


where,  consequently, 


VHo  =  2  S  {2[(kri)~3  +  jK)'2]  cos2  e. 


where 


-  [(kr.)“3  +  j(kr.)'2  -  (kr.)"1]  sin2  8. 

-jkr 

•  Je  (32) 

-r  —  jlCT* 

)  3  +  3j(kr.)  2  -  (kr.)-1j  e  1  sin  0. 

cos  0^  cos  0^  (33) 

r2  =  x2  +  d2  +  (z  +  h)2 

(34) 

sin  81  =  (x2  +  d2)l/2/ri 

(35) 

cos  0  =  ( z  +  h)/^ 

(36) 

cos  0X  =  x/(d2  +  x2)1/2 

(37) 

r2  =  x2  +  d2  +  (h  -  z)2 
sin  e2  =  (x2  +  d2j1/2/r2 
cos  «2  «  (z  -  h)/r2  . 


14 


(38) 

(39) 

(40) 
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Two  other  orientations,  orthogonal  to  the  first,  are  possible  and  may  be 
of  interest.  These  and  their  corresponding  image  theory  solutions  are  shown  in 
figures  8  and  9.  Expressions  corresponding  to  equations  (32)  through  (40)  are 
readily  obtained.  Thus,  the  results  for  an  arbitrarily  oriented  dipole  may  be 
obtained  through  suitable  superposition  of  these  cases.  If  the  dipole  in  figure  7 
were  oriented  at  angles  withrespect  to  the  cartesian  coordinate  system, 


m  = 

m 

cos 

0 

(41) 

z 

m 

m  = 

m 

sin 

e 

sin 

<t> 

(42) 

y 

m 

m 

m  = 

X 

m 

sin 

e 

m 

cos 

<t> 

m 

(43) 

is  obtained,  where  mz,  my,  and  mx  are  the  dipole  moments  to  be  inserted  in  the 
expressions  corresponding  to  the  cases  shown  in  figures  7,  8,  and  9,  respectively. 

It  should  be  emphasized  that  the  solution  obtained  holds  only  for  flu  region 
y  >  0, z  <  0  (the  ocean)  and  cannot  be  extended  successfully  beyond  that  region. 
Thus,  the  surface  current  induced  on  the  metallic  plane  can  only  be  obtained  for 
the  submerged  part.  As  indicated  earlier  in  this  section,  little  energy  pene¬ 
trates  into  the  air  region  and,  therefore,  the  induced  currents  above  the  line 
z  =  0, y  =  0  will  have  little  or  no  effect  on  the  final  result.  The  only  exception 
to  this  would  occur  if  the  dipole  were  located  close  to  the  coordinate  origin, 
where  it  would  be  in  the  immediate  proximity  of  the  conducting  plane  and  only 
just  under  the  surface.  In  this  case  substantial  errors  might  well  result. 


THE  INFLUENCE  OF  THE  AIR-OCEAN  INTERFACE  ON  A 
DETECTOR  ORIENTED  ORTHOGONAL  TO,  AND  AT  A 
DISTANCE  b  FROM,  THE  SOURCE 

The  results  of  the  previous  section  were  extended  to  allow  an  examination 
of  changes  in  signal  detection  due  to  the  presence  of  the  air-ocean  interface 
when  the  source-detector  system  approached  the  infinite  conducting  plate.  To 
this  end,  equations  (25)  through  (27)  become  pertinent  for  the  case  correspond¬ 
ing  to  figure  7  where  the  dipole  source  is  located  at  the  point  x  =  0,y  -  d, 
z  =  -h,  and  the  detector  is  located  at  the  point  x  =  b  cos  y  =  d  +  b  sin  $ jj, 
z  =  -h;  i.e.,  the  detector  is  presumed  to  be  located  in  the  same  horizontal 
plane  as  the  source  and  separated  from  it  by  a  distance  b.  The  direction  of  b 
makes  an  angle  of  with  the  y  =  0  plane. 
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The  expressions  for  (rj,0j,4>j)  become,  for  a  general  detector  location 

(xd,yd,zd). 


2  2 

r  =  x  + 
1  d 

(yd'd)2+(Zd  +  h)2 

(44) 

2  2 
r2  =  Xd  + 

(yd  ■  d)2  +  (Zd  ‘  hf 

(45) 

2  2 
r3  =  Xd  + 

(yd  +  d)2  +  (zd  -  h)2 

(46) 

2  2 
r4  =  Xd  + 

(yd+d)2+(Zd+h)2 

(47) 

sin  91=  [xd  +  (yd-d 

f]l/2/ri'  cos  ff1=(zd+h)/r1 

(48) 

8ine2  =  K+(yd-d 

)2]l/2/r 2’  C0S  ®2  =  (Zd“h)/r2 

(49) 

sin  e3=[xd  +  (yd  +  d 

)2]1/2/r3>  cos  e3  =  (zd-h)/r3 

(50) 

sin  ®4  =  [xd  +  (yd+d 

)T/2/r4’  cos  ®4  =  (Zd  +  h)/T4 

(51) 

sin  =  (yd  -  d)/[x2  +  (yd  -  d)2]1/2,  cos  ^  =  xd/[x2  +  (yd  -  d)2]l/2  (52) 

*2  =  \  (53) 

sin  «3  =  (yd  +  d)/[x2  +  (yd  +  d)2]1/2,  cos  <t>3  =  xd/[x2  +  (yd  +  d)2]l/2  (54) 

*4=*3.  (55) 


~™'  .  ■  1  \  - - 


WS«,pi.lBIJ,«vJ 
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If,  as  previously  mentioned,  the  detector  location  is  constrained  to  the 
horizontal  plane  contain. ng  the  source,  then  the  component  of  the  magnetic  field 
intensity  which  is  of  particular  interest  is  the  one  lying  in  the  direction  of  the 
baselength  b, 


H  =  H  cos  0,  +  H  sin  0,  , 
b  x  by  b 

since  in  the  absence  of  the  metallic  plane  and  the  air-ocean  interface  this  com¬ 
ponent  will  vanish.  The  vertical  component  at  the  detector  is 

H  =  H 
v  z 

and  the  component  normal  to  both  and  HQ  is 

H  =  H  cos  0  -  H  sin  <t>,  . 
n  y  b  x  b 

Figures  10,  11,  and  12  give  the  results  of  some  of  the  calculations  per¬ 
formed  on  this  system.  The  ordinates  represent  the  EMF  induced  in  a  coil  of 
radius  a  located  in  a  plane  normal  to  the  baseline  b  at  the  detector  location. 
The  EMF  is  normalized  with  respect  to  the  EMF  induced  in  the  same  coil  when 
located  with  its  axis  coincident  with  the  source  in  a  plane  normal  to  the  mag¬ 
netic  dipole,  and  is  expressed  in  dB's.  All  linear  quantities  are  normalized 
with  respect  to  the  electromagnetic  skin  depth  in  the  ocean.  The  variation  in 
the  signal  level  is  shown  as  the  dipok  source  approaches  the  conducting  plane 
at  various  constant  depths  from  0. 1  to  1.5  skin  depths.  The  effect  of  the  de¬ 
tector  orientation  with  respect  to  the  conducting  plane  is  shown  in  some  sample 
cases. 

It  is  evident  that  the  effect  of  the  air-ocean  interface  is  strong  and  masks 
the  effect  of  the  conducting  plane  at  shallow  depths  until  the  dipole  has  approached 
to  within  one-half  skin  depth.  At  greater  depths,  one  skin  depth  or  more,  the 
infinite  conducting  plane  is  first  "seen"  at  about  two  skin  depths.  The  three 
figures  correspond  to  baselengths  of  0.075,  0,150,  and  0.250,  successively, 
and  indicate  that  decreasing  the  baselength  increases  the  rate  at  which  the  sig¬ 
nal  changes  once  the  conducting  plane  has  been  "seen." 
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Figure  10.  EMF  Induced  in  Orthogonal  Detector  as  a  Function  of  Depth  and  Stand-off:  Baselength  0.  075  6 
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Figure  11.  EMF  Induced  in  Orthogonal  Detector  as  a  Function  of  Depth  and  Stand-off:  Baselength  0. 150  5 


Figure  12.  EMF  Induced  in  Orthogonal  Detector  as  a  Function  of  Depth  and  Stand-off:  Baselength  0.250  5 
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SUMMARY 

This  report  consists  oi  three  parts.  The  surface  current  distribution 
induced  by  a  magnetic  dipole  on  an  infinite  conducting  plane  was  investigated 
when  both  were  immersed  in  a  lossy  medium,  the  ocean.  It  was  found  that  the 
surface  current  distribution  did  not  extend  materially  beyond  a  circular  area  of 
about  one  skin  depth  in  radius  centered  on  the  point  on  the  conducting  plane 
nearest  the  dipole  source.  This  result  was  then  generalized  to  include  the  case 
where  the  effect  of  an  air-ocean  interface  was  included.  In  the  third  and  final 
part,  the  signal  induced  in  a  detector  coil  was  obtained  for  a  vertical  magnetic 
dipole  in  the  ocean  in  the  presence  of  an  infinite  vertical  perfectly  conducting 
plane  and  a  horizontal  air-ocean  interface.  The  spacing  between  the  detector 
and  source  was  varied  in  the  same  horizontal  plane,  while  the  detector  coil 
was  oriented  orthogonally  to  the  magnetic  dipole  and  its  axis  was  in  line  with 
the  baselength.  It  was  found  that  the  air-ocean  interface  affected  the  signal 
severely  and  masked  the  presence  of  the  conducting  plane  for  shallow-depth 
runs.  At  greater  depths,  as  expected,  the  air-ocean  interface  effect  is  rapidly 
lost.  Appendixes  A  and  B  contain  the  computer  programs  used  for  computations 
pertaining  to  the  first  and  third  parts  of  this  report. 
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Appendix  A 

CURRENT  DISTRIBUTION  ON  AN  INFINITE  PLATE  AS  A 
FUNCTION  OF  MAGNETIC  DIPOLE  DEPTH  FROM  THE  PLATE 


C  CURRENT  DISTRIBUTION  J  SUB  PHI  ON  AN  INFINITE  PLATE  AS  A 
C  FUNCTION  OF  MAGNETIC  DIPOLE  DEPTH  FROM  PLATE.  S  EQUALS 
C  S  OVER  DELTA,  D  EQUALS  D  OVER  DELTA  WHERE  DELTA  IS  THE 
C  SKIN  DEPTH,  D  IS  A  MEASURE  OF  THE  DIPOLE  DEPTH,  S  IS  A  MEAS- 

C  URE  OF  THE  DISTANCE  FROM  THE  POINT  ON  THE  PLATE  DIRECTLY 

C  ABOVE  THE  DIPOLE. 

READ  1 ,  DMIN,  DMAX,  DD,  SMAX,  DS 

1  FORMAT (5  FI  0.4) 

D-DMIN 

3  S=0. 

DENOM=(((l.+(2.**.5))**2+((4.-(2.**.5))**2))**.5 

2  RATIO=S*D/((S*S+D*D))**2. 5) 

SUM=(S*S+D*D)**.5 

TOP=((l .  +SUM)**2+(2.  *SUM*SUM-SUM)**2)**.  5 

EX=-SUM 

FF=EXPF(EX) 

TOTAL=4.  *(2.  **.5)*RATIO*FF*TOP/DENOM 
PRINT  4,  D,  S,  TOTAL 

4  FORMAT(5H  D=  ,E10,4,5H  S=  .E10.4.14H  J  SUB  PHI  =  ,E15.8) 

S=S+DS 

IF(S-SMAX)2,2,  5 

5  D=D+DD 
IF(D-DMAX)3, 3,  6 

6  END 


Note:  The  current  distribution  has  been  normalized  such  that  the  output  is  in 
the  form  J^(d,s)/J^(a,  j). 
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Appendix  B 


SIGNAL  INDUCED  IN  A  DETECTOR  COIL  AS  A  FUNCTION  OF  THE 
SUBMERGED  DISTANCE  FROM  AN  INFINITE  CONDUCTING  VERTICAL  PLATE 


C  SIGNAL  INDUCED  IN  DETECTOR  COIL  AS  A  FUNCTION  OF  SUB- 
C  MERGED  DEPTH.  DISTANCE  FROM  AN  INFINITE  CONDUCTING 
C  VERTICAL  PLATE.  BASE  LENGTH  AND  ANGLE  OF  BASE  LENGTH 

C  TO  VERTICAL  FROM  THE  INFINITE  PLATE  TO  THE  DIPOLE  ELE- 

C  MENT.  DIPOLE  ORIENTATION  IS  VERTICAL  AND  THE  NORMAL  TO 
C  THE  DETECTOR  COIL  LIES  IN  THE  HORIZONTAL  PLANE.  FURTHER 

C  PARAMETERS  ARE  -A  THE  COIL  RADIUS,  N1  NUMBER  OF  TURNS 

C  ON  PRIMARY  COIL,  N2  NUMBER  OF  TURNS  ON  THE  DETECTOR 

C  COIL.  ALL  LINEAR  DIMENSIONS  ARE  NORMALIZED  W.R.T.  THE 

C  SKIN  DEPTH  IN  SEA  WATER.  PHI  IS  TO  BE  ENTERED  IN  DEGREES. 
REAL  N1,N2 

DIMENSION  HH(10),  BB(10),  BPHI(IO) 

HH(1)  =  .  1 
HII(2)  =  .25 
HH(3)  =  .  50 
HH(4)  =  1. 00 
HH(5)  =  1. 50 
BB(1)  =  .  075 
BB(2)  =  .  150 
BB(3)  =  .250 
BPHI(1)=0. 

BPHI(2)=30. 

BPHI(3)=60. 

BPHI(4)=90. 

DMAX=2.000 
A =.025 
Ni=i. 

N2  =  l. 

DO  6  1=1,  5,1 
H=HH(I) 

DO  8  J=l, 3, 1 
B=BB(J) 

DO  7  K=l,  4, 1 
BHI=BPHI(K) 

PRINT  2,DMAX,H,B,BHI,A,N1,N2 
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2  FORMAT(7H  DMAX=  E10.3.4H  H=  E10.3.4H  B=  E10.3.6H 
1PHI=  E10.3.4H  A=  E10.3.5H  Nl=  E10.3.5H  N2  =  E10.3) 

PRINT  3 

3  FORMAT (40H  D  EMF  EMFDB) 

BHI  =  (BHI/180.)*3. 1415927 

DO  9  M=l, 40 
D  =  .05*(M-1) 

R3  =  (4.  *H*H+B*B)**.  5 

R4=(R3*R3+4.  *D*D44.  *D*B*COSF(BHI))**.  5 

CS3=2.*H/R3 

SI3=B/R3 

CS4=2.*H/R4 

SI4=(4.  *D*D+B*B+4.+D*B*COSF(BHI))*'\  5/R4 
BHI4=ATANF(B*SINF<BHI)/(2.*D+B*COSF(BHI))) 

DENOM  =  (-A*COSF(A)+{l.+A)*SINF(A))*EXPF(-A)/(A*A) 
T3=A/(R3*R3*R3) 

T4=A/(R4*R4*R4) 

S3=(l .  +R3)*SINF(R3) 

S4=(l.  +R4)*SINF(R4) 

TOPl=T3*(-R3*COSF(R3)+S3)*EXPF(-R3)*CS3 

TOP2 = (T3/2 .  )*  (-R3*(l .  +2 .  *R3)*COS  F  (R3)+S3)  *EXPF  (-R3)*SI3 

TOP3=T4*(-R4*COSF(R4)+S4)*EXPF(-R4)*CS4 

TOP4=(T4/2 .  )*(-R4*(l.  +2.  *R4)*COSF(R4)+S4)*EXPF(~R4)*SI4 

CS5=COSF(BHI-BHI4) 

TOP=TOP1*S13+TOP2*CS3+TOP3*SI4*CS5+TOP4*CS4*CS5 

EMF=TOP/DENOM 

EMF=ABSF(EMF) 

EMFDB=20.*ALOG(EMF)/2. 302  585 
PRINT  4,  D,  EMF,  EMFDB 
4  FORMAT(3E15, 5) 

9  CONTINUE 

7  CONTINUE 

8  CONTINUE 
6  CONTINUE 

END 
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